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^ ■ Abstract 

■ We analyze the transverse momentum distribution of J/ ijj mesons produced in Au + Au collisions 
at the top RHIC energy within a blast-wave model that accounts for a possible inhomogeneity 

^ \ of the charmonium distribution and/or flow fluctuations. The results imply that the transverse 

momentum spectra of J/V', (p ^"^^ ^ hadrons measured at the RHIC can be described well if kinetic 
freeze-out takes place just after chemical freeze-out for these particles. 
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I. INTRODUCTION 



Charmonium production is considered as important messenger of the deconfinement tran- 
sition in heavy-ion colhsions. In the original investigation the suppression of charmonium 
was proposed as a signal of the quark gluon-plasma (QGP) that caused the melting of char- 
monium [l| (for very recent reviews see Refs. 0, 3]). Alternatively, the charmonium yield is 
described as been caused by statistical hadronization at the phase boundary [4] E In the cor- 
responding statistical hadronization model (SHM), the J /ip multiplicities at SPS and RHIC 
energies are well described for different centrality and rapidity bins; the basic points of the 
SHM and its current development are presented in reviews and references therein. 

The SHM assumes that the charm quarks are produced in primary hard collisions and that 
their total number stays constant until hadronization. A crucial hypothesis of the model is 
thermal equilibration of charm quarks in the QGP, at least near the critical temperature T^. 
The charmonia are then all produced near the phase boundary!^ 

Because the momentum spectrum of J / ip is nearly frozen out at T^. (the typical cross sec- 
tion of J /if) with CO moving hadrons is at most a few millibarns |9;]), the measured momentum 
spectra of J /ip mesons contain valuable information on details of the hadronization process 
in the strongly interacting medium. For instance, if formation of J /ip proceeds as assumed 
in the SHM, then the momentum spectrum of J/ip is expected to be locally equilibrated 
in the vicinity of T^. On the contrary, if the superdense matter formed in A+A reactions 
contained a significant portion of J/ip mesons produced at the early stage of the process in 
hard nucleon-nucleon collisions, then the final J/ip spectra could exhibit two components, 
as, for example, in the model [lO|: direct J/ip and secondary thermal J/ip produced at the 
stage of recombination of charm quarks. The shape of J /ip momentum spectra will then be, 
in fact, nonthermal. Also, in the latter approach, the multiplicities of hadrons with open 
and hidden charm do not coincide, generally, with those predicted in the SHM. 

Recently the PHENIX Collaboration published results for J/ip production versus cen- 
trality, transverse momentum and rapidity in A+A collisions at ^/snn = 200 GeV jll l. 
The analysis of these data in Ref. ^| by means of the hydroinspired blast-wave model 12 1 
^ Possible thermal features in charmonium and, in particular, J/?/; production in Pb+Pb collisions at SPS 
energy were first pointed out in Ref. 5[. However, thermal production of charm quarks is negligible at 

RHIC energy and below because of the large charm quark mass [4 1 . 
^ Perhaps slightly above T^. the upper limit for the dissociation temperature was recently estimated as 

1.2Tc 



demonstrated that a simple blast-wave model parametrization of the freeze-out does not 
fully account for the results measured by the PHENIX Collaboration on J / ip transverse mo- 
mentum spectra despite the success of the SHM in the description of multiplicities. On the 
contrary, the fit based on the two-component model [l^ with the blast-wave parametrization 
of transverse momentum spectra for secondary J / ip describes well the momentum spectra of 
J/ip isj]. However, the portion of primordial J/ip that were not melted in the QGP seems 
to be too big in this model: about 50%. Note that a good description of J/ip spectra can 
also be reached without the assumption of thermal equilibration of the secondary charmo- 
nia. If recombination of charm quarks is based on the coalescence mechanism (for review 
see, e.g., Ref. and references therein), then the secondary J/ip are not thermal even if 
charm quarks are considered to be thermalized ones [15|] at the hadronization stage. The 
latter scenario is realized, for example, in Ref. [l^ where the recombination coalescence 
component of J/ip spectra was calculated based on the kinetic coalescence model of Re 
17l | . However, although the recombination contribution matches the low pt region well 
significant primordial J/ production is again required to describe the higher pt region 
In the present article we advocate the local equilibrium picture of J/ip production at 



s 



a. 



m 



the freeze-out stage near in accordance with the statistical hadronization model 

I?]. While multiplicities in this model were described in Refs. 0, 0], our aim here is to 
demonstrate that transverse momentum spectrum slopes of J/ip mesons at midrapidity 
can be described by a hydroinspired parametrization without a significant primordial J/ijj 
contribution: the latter is questionable taking into consideration recent results Q| as for J/ijj 
melting in the QGP. The article is organized as follows. In Sec. II, we present the results 
of blast-wave model fits to describe transverse momentum spectra of and Q hadrons at 
RHIC and discuss the problems with the fit of J/ip spectra. In Sec. HI, we briefiy analyze 
the main assumptions of the blast-wave model, to clear up the reasons for the shortcomings 
in the blast-wave model description of J ftp transverse momentum spectra, and show that a 
consistent description of the (p, Q, and J ftp transverse momentum spectra can be obtained 
within a hydroinspired parametrization of transverse momentum spectra that accounts for a 
possible inhomogeneity of charmonium distribution and/or fiow fiuctuations. We conclude 
in Sec. IV. 



See also Ref. 



18| , where the J/ip transverse momentum distribution was calculated in a transport model 



with both initial production and continuous regeneration of charmonia. 
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II. BLAST- WAVE MODEL FITS TO TRANSVERSE MOMENTUM SPECTRA 
OF (p AND n AT RHIC AND TRANSVERSE MOMENTUM SPECTRA OF J/ip 



Let us start with the assumption that the momentum spectra of J/ip, as well as and 
Q hadrons, are frozen near the phase boundary in RHIC-energy collisions at T 
This value is close to the temperature of chemical freeze-out in these collisions 
also note that Tc as obtained for zero-net-baryon-density matter from the most recent lattice 



60 MeV. 
3- We 



QCD calculations is close to this value 



2l| . We choose (p and Q particles for comparison 



with J/ip mesons because these particles are thought to have very small hadronic cross 
sections and their spectra are not expected to be significantly distorted by feed-down from 
resonance decays in relativistic heavy-ion collisions. The assumption of coincidence of the 
kinetic freeze-out of these particles with their chemical freeze-out was supported in Ref. 22 1 
where it was demonstrated that the blast-wave model fit to the transverse mass spectra of 
these particles at RHIC yields T ^ 160 MeV along with an average fiow velocity (v) ~ 0.45c. 
The latter is lower than that for hadrons strongly interacting in the hadron resonance gas, 
and, thereby their chemical freeze-out precedes kinetic freeze-out Coincidence of thermal 
and chemical freeze-out also underlies the hydroinspired model presented in Ref. 24 1. 

It is noteworthy that the inverse slope parameters for Q and transverse momentum 
spectra do not vary significantly with collision centrality 25l-l27l|. This weak centrality 



dependence can be explained by the independence of the chemical freeze-out temperature 
from the collision centrality that was observed in the chemical equilibrium model fit to 



particle number yields [20|, 



281]. This experimental observation also supports the coincidence 



of kinetic and chemical freeze-out for these particle species. Note that the insensitivity of the 
chemical freeze-out temperature to centrality is natural 29| if chemical freeze-out happens 
at the phase boundary [3(| (see also the recent review article 31]). 

To calculate transverse momentum spectra of J/ip, 4> and Q particles in central collisions 
we utilize the blast-wave formula 121 ] . 



dN 



Prdprdy 



oc rriT 



R 



drrlo 







Pt sinh yT{r] 



T 



tut cosh 1/7- (r) 
T 



(1) 



^ A similar conclusion was also obtained in Ref. [23] for SPS 158^ GeV energy collisions where blast-wave 
fits to J/ip and ip' mesons and 51 hyperon spectra in Pb+Pb collisions were performed to show that a good 
description of these particle spectra by the blast-wave formula with the temperature associated with the 
chemical freeze-out and a relatively low transverse collective velocity can be obtained. 



which follows from boost-invariant parametrizations of freeze-out at r = y/t'^ — z'^ = const 
under the assumption of homogeneous particle number density until r = R, where R is 
the transverse size of the system. Here Iq and Ki represent modified Bessel functions, and 



rriT = \/Pt + We utilize a constant temperature T across the freeze-out hypersurface 
and a linear transverse rapidity profile, t/T = y^^^r/R, where y^^^ = yriR) is the maximum 
transverse rapidity. The collective radial expansion velocity v is then v = tanhy^-, and 
"Wmax = tanh y^^ is the velocity at the surface. 

Because we are interested in the description of the transverse spectrum slope, it is con- 
venient to substitute r/R with x in Eq. ([1]) to obtain 

dN , ^ f pTsmh.(ax)\ /m-r coshfax) \ 

(xtut xdxio — — — ^ Ki — — — ^ — ^ , (2) 



PTdprdy Jq \ T 

with a = y^^"^. 

We check first that assuming equal chemical and kinetic freeze-out temperatures is a 
good approximation for and (^7 = Q'^ + Q^) particle transverse momentum spectra 
at midrapidity. The corresponding experimental data for Au+Au collisions at a/s/vtv = 
200 GeV for Q baryons in the 0% — 10% centrality bin are taken from Re: 
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|25| (STAR 
(PHENIX 



Collaboration), and those for mesons at 0% — 10% centrality from Ref. 

Collaboration). As shown in Fig. [H the shapes of the fl and spectra are reproduced well 

by the blast-wave formula with fairly realistic parameters at chemical freeze-out: T = 160 

MeV, a = 0.75, implying v^^^ ^ O.GSsEI We also show in Fig. [1] the exponential fit to 

momentum spectra of mesons in the 0% — 10% centrality bin by the STAR Collaboration 

271], which we normalized to PHENIX data points. The STAR Collaboration demonstrated 

that 0-meson momentum spectra are fitted well by an exponential function, 

d'^N 

exp (-mr/Texp) (3) 



prdprdy 

where the slope parameter Texp ~ 359 MeV for the 0% — 10% centrality bin. Note that 
the dN/dy results for mesons extracted by the PHENIX and STAR Collaborations in 
this centrality bin are not in agreement, while, as explicitly shown in Fig. [H there is no 
discrepancy between the measured spectrum slopes. 

^ Note that relatively high collective transverse velocities for T — 160 MeV can appear in the course of 



the evolution |32| as a result of the development of initial transverse velocities at the prethermal partonic 
stage (with subsequent rapid thermalization) and viscous effects. The latter leads to a decrease in the 
longitudinal velocity and increase in transverse velocity. 
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FIG. 1. Transverse momentum data (in arbitrary units) at midrapidity of (j) mesons and ($7 + 
0,^) hyperons in central Au + Au collisions at y^J/vJV = 200 GeV RHIC energy along with the 
results of the blast- wave model, Eq. Q, fit which are indicated by the solid lines, the dashed 
line represents the exponential, Eq. ([3|), fit to momentum spectra of mesons by the STAR 



Collaboration [27 
Q are from Refs 



which we normalized to the PHENIX data points. Experimental data for (/) and 



261] and 



25[ respectively. Statistical error bars are shown for (p mesons, and total 



(statistical and systematic) error bars are presented for Q hyperons. 

The results for momentum spectra of J/ip with the same set of parameters are shown 
in Fig. [21 together with the data measured by the PHENIX Collaboration ll| for two 



centrality binsj^ One notices significant deviations for the data points at low momentum, 
leading to a excessively high effective temperature. 

To see what blast-wave model parameters are required to describe the low-p^ region of 
J/ip spectra, we performed a blast-wave model fit to the low- momentum part of J/ip spectra 
at the same temperature, T = 160 MeV, but considered the maximum transverse rapidity 
a in Eq. ([2|) to be a free parameter. The results are demonstrated in Fig. [21 where one can 
see that the blast-wave model with a = 0.55 (implying t>max ~ 0.5) yields a good description 



We show the data points of noncentral collisions (20% — 40% centrality bin) just for comparison. 
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FIG. 2. Transverse momentum data for J/tp mesons (in arbitrary units) at midrapidity in Au + Au 
collisions at ^/sNN = 200 GeV RHIG energy along with the results of the blast-wave model, Eq. 
([2]), fit with the same parameters as in Fig. [1] (solid lines). Dashed lines represent the blast- 
wave model, Eq. ([2|), fit results with the "soft" transverse flow parameter: see text for details. 
Experimental data for J/-!/; are from Ref. ll[ for two centrality bins. Statistical and systematic 



uncorrelated errors are represented as error bars. 

of the low-momentum part of the J/ip spectrum. These findings, at first sight, might imply 
a contradiction with the blast-wave model fits to (p and Q transverse momentum spectra 
where the same temperature but another surface velocity was obtained. Consequently, the 
question arises: Can J/tp, (j) and VL transverse momentum spectra be coherently described 
within a hydroinspired parametrization of freeze-out with reasonable parameters? 



III. HYDROINSPIRED PARAMETRIZATION OF J/ip TRANSVERSE MOMEN- 
TUM SPECTRUM FREEZE-OUT AT THE RHIC 

To understand why a blast-wave model, Eq. ([1]), description with parameters fitted to 
and VL transverse momentum spectra does not reproduce the J /ip spectra properly, let us 



look carefully over the main assumptions used in the comparison to spectra by employing 
Eq. ([1]). They are the following: the freeze-out isotherm is r = const and the transverse 
rapidity profile at the isotherm is linear. Further assumptions are on-mass shell distribution 
functions, absence of flow fluctuations, a homogeneous particle number density, n = y/nPm^, 
at freeze-out, and absence of resonance feed-down. Here we briefly discuss these assumptions 
to clarify whether giving up debatable ones can result in consistent description of 0, Q, and 
J/ip transverse momentum spectra if kinetic freeze-out takes place just after chemical freeze- 
out for these particles. 



First note that a hypersurface of constant proper Bjorken time, r = y/W--z^ = const 
in general does not correspond to any isotherm, and a linear transverse rapidity profile is 
questionable. However, based on estimates of corresponding quantities in viable scenar- 
ios of hydrodynamic evolution and results of correspondingly modified blast-wave model 



calculations (see, e.g., Ref. [20| for blast-wave model calculations with different velocity 
profiles), one can conclude that the spectra discussed are little affected by the shape of 
the flow velocity profile and by deviations of the chemical freeze-out isotherm from the 
r = const hypersurface. We hence keep these assumptions unchanged. Actually, the cor- 
responding hypersurface and velocity profile depend on the whole dynamical history of the 
hydrodynamically expanding system - initial conditions, equation of state, etc. - and can be 
calculated in a true dynamical model only. It will be the subject of follow-up work. 

Another assumption, namely, the on-mass shell approximation of the distribution func- 
tion, may be questioned for finite-width particles such as J/ip, although the free width of 
the J/ip is less than 100 keV. The main reason for the higher effective temperature of J/ijj 
spectra in the low-momentum region (see Fig. |2]) compared with the effective temperatures 
of and Q spectra in the same region (see Fig. [1]) is the larger mass of the J/ip particle. An 
explanation of this fact can be based on the simple approximation of the integrand in Eq. 
(|2]) . Then one may surmise that an appropriate off- mass-shell hydroinspired parametrization 
of J/ ip spectra that includes low invariant mass contributions could improve the description 
of the low-momentum part of the spectra. Most naturally, such a parametrization - namely, 
a Breit-Wigner formula for the spectral function instead of the S function approximation - 
appears because of the finite width of J/ip meson. However, the corresponding calculations 
of the J/ip spectra with the same hydrodynamical parameters of the blast- wave model as 
for Q and particles do not result in a significant improvement in the description of J/ijj 
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spectra if experimental conditions are properly taken into account. Namely, as reported by 



the PHENIX Collaboration ll|, the J /if) spectra at midrapidity are measured by counting 
lepton pairs in the invariant mass range 2.9 GeV < M <3.3 GeV, and such a narrow invari- 
ant mass window accompanied by the rather low decay width, Ftot = 93.2 x 10^^ GeV, does 
not result in noticeable deviation from on-mass shell calculations. 

The next assumption, namely, a neglect of flow fluctuations in the blast-wave model, 
Eq. (Il]), is also debatable. Owing to the finite size of the systems, large fluctuations are 
expected in the initial stage of the nuclear collisions, even for a fixed impact parameter, and 
various event generators do show such effects. While a very early thermalization in rela- 
tivistic heavy-ion collisions is doubtful now 33|] (for review see also Ref. and references 



therein), the initial inhomogeneities can result in fluctuating initial conditions for the sub- 
sequent hydrodynamical expansion. Because the hydrodynamical equations are nonlinear, 
the event average of any hydrodynamical parameter (in other words, average over solutions 
of hydrodynamical equations) is quite different from that for a smooth initial configuration 
and results in large differences in spectra from hydrodynamical calculations with averaged 
initial conditions. This fact was pointed out and explicitly demonstrated in event-by-even ; 
hydrodynamic simulations based on smoothed particle hydrodynamics code (see Ref. jsSj 



for review and Ref. 



for recent results). Fluctuations of transverse flows in noncental 



collisions also are considered in Ref. 



37|. 



Therefore, the collective flow fluctuations are rather natural and can have some effect 
on the particle momentum spectra. To study it, let us consider a blast- wave model with 
bulk matter flow fluctuations. The transverse spectra of the blast-wave model averaged over 
ensemble of the fluctuations are assumed to have the form: 

cx niT / daG{a) / dxxlo — ^ Ki — ^ — . (4) 



Prdprdy J^^.^ Jq \ T J \ T 

We consider here two simple and, in some sense, opposite cases of the distribution of a. 
First, we assume a distribution that is flat in a, that is, G{a) = 1, with the lower and 
upper limits amin and ctmax considered to be free fit parameters. Second, we assume a 
Gaussian form for the distribution of hydrodynamical velocities. With f max = tanh a, we 
obtain: G(a) = exp (tanha-tanh«°)^ j ^ (^- ("— with «min = (v^^ = 0) and 
«max = oo (fmS = ^) ■ Notc that, for such a distribution, there is no explicit cutoff of high 
a contributions; the convergence of the integral takes place because the main contribution 



at high a happens at small x ~ 1/a, leading to a cutoff factor l/a^. The results for 
G{a) = 1 with a^in = 0.23 (^^"^ ^ 0.23) and a^ax = l-O ^ 0.76), and for G{a) = 

exp ^— ^^'^"'^""^2'^"^" ^ j with a° = 0.55 (fmax ~ 0-5) and 6 = 0.18 are displayed in Figs. 
[3] and m for Q, 0, and J/ip transverse momentum spectra, respectively. One can see from 
these figures that the blast-wave model supplemented with appropriate fluctuations of the 
hydrodynamical flow yields a significantly improved description of the shapes of transverse 
momentum distributions for (p, Q, and J/ip particles. 
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FIG. 3. The same as Fig. [H but with the blast- wave model fit accounting for flow fluctuations, 
Eq. The corresponding results are indicated by solid lines (flat distribution in hydrodynamical 
velocities) and dashed lines (Gaussian distribution in hydrodynamical velocities). 



Another assumption used in Eq. ([T]), namely, that of homogeneous particle number 
density across the system at the isotherm hypersurface, is rather natural for bulk-matter 
particles consisting of quarks that are produced in QGP, as well as in the initial collisions, 
but can be questioned for heavy rare particles like J/ip. The reasons are the following. 
Because the charm quarks are produced in primary binary nucleon-nucleon collisions only 
and there is no charm quark production during the evolution of the system, the initial 
distribution of charm quarks follows the distribution of initial binary collisions. Then the 
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FIG. 4. The same as Fig. [2] but with the blast-wave model fit accounting for flow fluctuations, 
Eq. The corresponding results are indicated by solid lines (flat distribution in hydro dynamical 
velocities) and dashed lines (Gaussian distribution in hydro dynamical velocities). 

charm quark density can be more strongly peaked in the center of the system than the density 
of light-flavor and strange quarks that initially follow the participant density distribution!^ 
Hence the effective transverse size of a volume occupied by J/ip mesons at the hadronization 
hypersurface can be less than that for bulk-matter particles. 

For simplicity, we consider the particle number density of J/ip in two extreme scenarios: 
first, the J/ip mesons are assumed to follow the density distribution of bulk- matter parti- 
cles, as in fact considered in the previous section, and second, we assume that the spatial 
extension of J/ ip in the transverse direction coincides with the initial spatial distribution of 
charm quarks determined by the geometry of binary nucleon-nucleon collisions. The latter 
is examined here. 

The initial spatial distribution of charm quarks in the transverse direction as determined 
by the transverse distribution of binary nucleon-nucleon collisions in Au+Au collisions, and 

^ It was recently demonstrated [ssj that a difference in the initial distribution of charm quarks compared 
to the light flavors survives the hydrodynamical evolution and is still significant at decoupling, thereby 
resulting in the reduction of the mean transverse momentum of the D-meson spectra. 
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at impact parameter 6 = 0, which we assume here for simphcity, is 

n{r) oc Tj(r), (5) 
where the nuclear thickness distribution Tyi(r) is obtained from a Woods-Saxon distribution, 



TA(r) 



+00 



dzp ^Vr^ + . (6) 

The Woods-Saxon distribution of nucleons in the Au nucleus is defined by the formula 

v^r2~+^ — c 



with c = 6.38 fm, a = 0.535 fm 



1 + exp 



(7) 



and po given by the normalization condition. 

The transverse momentum spectrum of the blast-wave model of J/ tp distributed according 

to this initial spatial distribution of charm quarks is 

dN , ^9, ^ f PTsmh(ax)\ f rriT cosh(ax)\ , , 

oc ,nr I .i.Tl(.R)k [?^-^ ) ( ^ j . (8) 

Then, before calculating the transverse momentum spectra of J/ip by means of Eq. (IHl), 
one needs to define the transverse size of the bulk-matter distribution R. We estimate this 
value to be between 7 and 8 fm, and Fig. [5] presents the J / ip spectrum calculated with 
the hydroinspired parametrization, Eq. ([8]), for both R values. A good description of the 
spectrum at low pt spectra is obtained, this is not surprising because, as noted at the end of 
the previous section, to fit the low-p-p part of the J /ip spectrum one needs to utilize, for J /ip 
mesons, a lower value of the maximum transverse velocity than for bulk-matter particles at 
the same isotherm. 

The last assumption of the blast-wave model, Eq. ([T]), discussed here is the neglect of 
the resonance feed-down. While the measured inclusive particle spectra in general contain 
contributions from resonance decays, for the spectra discussed in this article the resonance 
feed-down is rather small and one can expect that generalization of the blast-wave model to 
include the resonance feed-down does not seriously influence the fltted parameters!^ In any 
case, feed-down contributions to J /if) spectra will, except for LHC energies where feed-down 
from i?-meson decay is likely signiflcant, coming from excited charmonia. Although this 



contribution can be up to 40 % {2], we expect little modiflcation of the spectral shapes. 



See also Ref. 



20| where it was demonstrated that resonance decays have no significant effect on the kinetic 



freeze-out parameters extracted by means of the blast- wave model. 
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FIG. 5. The same as Fig. O but with the blast-wave model fit accounting for inhomogeneity of 
the J /if) distribution, Eq. (I8|), with i? = 7 f m (solid lines) and i? = 8 f m (dashed lines). 

as the charmonia masses are all very large compared to masses of hadrons made of light 
and strange quarks. A comprehensive treatment of resonance feed-down within a dynamical 
hydrokinetic model of the evolution of the fireball and of particle production will be provided 
in a forthcoming work. 

IV. RESULTS AND DISCUSSION 

We have analyzed central J/ip., (f), and Q transverse momentum spectra at midrapidity 
within a blast- wave model and demonstrated that, while the high-pT- part of the spectrum of 
J/ip is well described by fit parameters obtained for a description of (p and Q spectra, the low- 
Pt experimental transverse momentum spectrum of J/ip exhibits an effective temperature 
lower than that calculated with the fit parameters for and Q spectra. In our opinion, 
the reason for this is grounded in some specific assumptions of the blast-wave model rather 
than in significant nonthermal contributions to J/ip spectra. To demonstrate this we have 
developed and presented a generalized blast-wave model. 
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First, we have demonstrated that a blast-wave model modified to account for collective 
flow fiuctuations results in good agreement with data over a wide pt region for J/ip, (j) 
and Vt transverse momentum spectra. However, good agreement with J/ip data is reached 
only if rather high flow fluctuations take place near the chemical freeze-out isotherm. Note 
that in a very recent paper [40] it is estimated that much lower flow fluctuations, following 
from the fluctuations of the transverse size of the initial source, are consistent with recently 
observed Wm event-by-event fluctuations of the average transverse momentum of bulk- matter 

n 

particles in Au+Au collisions at RHIC energiesjj If this result is conflrmed, it rules out the 
interpretation that the low-p^ peculiarities of J / ip transverse momentum spectra are a result 
of hydrodynamic flow fluctuations. 

Another generalization of the blast-wave model considered in this article is based on the 
fact that, because charm quarks initially follow the density distribution of primary binary 
nucleon-nucleon collisions, the shape of J/ip density at freeze-out differs from the particle 
density of uniform bulk matter because of a higher concentration of charm quarks in the 
center of the system at the hadronization hypersurface. The spectrum of J/V", calculated 
with the density distribution of primary binary collisions and with a conservative estimate 
of the transverse size of the volume occupied by the bulk matter, demonstrates a fairly good 
agreement with 0% — 20% J /ip transverse momentum data until pt = 3.5 GeV, and the curve 
calculated is below the data points for higher px- Note that high-p^- J/tp can be emitted 
mostly by corona and can be, in fact, the result of charmonium production in nucleon- 
nucleon collisions (see Ref. Q]). Because the fraction of participating nucleons contained in 
the corona region increases when the centrality decreases, in noncentral collisions one can 
expect a higher effective temperature of the J/ip spectra than in central collisions; indeed 
Fig. [5] shows this effected 

Finally, we conclude that the present data at the top RHIC energy are compatible with 
the picture in which J/ip, 0, and Q momentum distributions are frozen simultaneously 
with chemical composition at T = 160 MeV. We have demonstrated that a consistent 
description of the 0, Q, and J/ip transverse momentum spectra can be obtained within 
^ It is worthy of note, also, that the magnitude of event-by-event elhptic flow fluctuations in Au+Au 



collisions at ^/snn = 200 GeV 



42| was found to be in agreement with predictions based on spatial 



fluctuations of the participating nucleons in the initial nuclear overlap region. 
^'^ The relatively high corona contribution is, perhaps, the reason why the _D'^-meson spectra measured by 

elatively high radial 



the STAR Collaboration in the 0% — 80% centrality bin 43] can be described by the ri 



flow velocity at the surface, 0.6, and the freeze-oi^ temperature, 170 MeV (see Ref. 



a hydroinspired parametrization of common freeze-out for these particles if nonhomogeneity 
of the charmonium distribution and/or significant flow fluctuations are assumed. 
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